Comment on 'Evolution of the unidirectional electromagnetic pulses in an anisotropic 

two-level medium' 



N.V. Ustino\0 

(Dated: August 28, 2008) 

Recently, Zabolotskii [Phys. Rev. E 77, 036603 (2008)] presented the Lax pair for a version of the 
reduced Maxwell-Bloch equations. This version was derived under considering the unidirectional 
propagation of a two-component electromagnetic pulse through an anisotropic two-level medium. 
We demonstrate that his derivation contains essential omission, which led to the wrong version of 
the reduced Maxwell-Bloch equations. We also point out that the Lax pair for correct version of 
these equations had been known in the most general anisotropic case. 
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PACS numbers: 41.20.-q, 02.30.1k, 42.50.Md, 42.65.Tg 

A version of the reduced Maxwell-Bloch (RMB) equa- 
tions was derived in Ref. [l| as follows. The Hamiltonian 
of the interaction of two-component electromagnetic field 
with anisotropic two-level medium was taken in the form 

H = ^fa,-{d^E' x + d^E' v )a lx 

-(d<£>E' x + d$E' y )a 22 - (d xx E' x + dvyE'^ax 

-(dy X E' x +dyyE' y )a 2 , (l) 

where &k (k — 1,2,3) are the Pauli matrices, an — 
diag(l,0), (T22 = diag(0, 1), E' x and E' y are the trans- 
verse components of the electric field, loq is the frequency 
of the transition, h is the Plank's constant. Equivalent 
representation of H is 



H 



-— a 3 - d x E x 



dyE' y , 



(2) 



where the matrices d x and d y corresponding to the op- 
erators of the dipole moment components are defined as 
given 
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Coefficients d z k ) (k = 1,2, s = x,y), d sp (s,p — x,y) of 
these matrices are assumed to be real constants. 

New variables were introduced by the transformation 



where 
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fix — \J dl x + d 2 xy , S y — yj d yx + d y 
Then, Hamiltonian |T]) reads as 



H 



Here 



hujQ 



-<T 3 



{p^E x +p^E y )a 11 



2 \rx — * ■ ry 

-{p {2) E x + p y 2) E y )a 22 - 6 x E x a 1 - 6 y E y a 2 . 



d (k) d -d (k) d 

Wzx Mini Uizil tX-j t 



Px 



(fe) _ c ^ u yy u yx (k) _ ? 

X X p , Py — 0., 



wWw -d [k) d 



(/c 1,2), Pq — d xx dy y d X ydy X . 

The dynamics of the two-level medium was described 
by the von Neumann equation on the density matrix p, 

ih^=[H,p], (4) 

with 



p= PU 912 ). 

\ P21 922 J 

Eq. ([4]) was rewritten in the terms of the components of 
the Bloch vector 

„ _ Pl2 + p2l „ _ P\2 — P21 „ _ Pll — P22 

bx ~ 2 ' y ~ 2i ' 2 ~ 2 
and the dimensionlcss variables 

_ dpE x _ doEy 

&x — ~~ Z ; &y 



frhjQ 



where 



d = y + + (pP - Px 2) y + ( 



(1) (2)\ 

Py -Py ) > 



as given 



dS. 
cV 



o , — (1 rrix£x m y£y)^y + fJ,y£yS Z , (5) 



dSy 

dr' 



{m x £ x + m y £y ~ \)S X - p x £ x S z , (6) 
dS 



0t 



f f^X Sy /i y Sy S x 



(7) 
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Here r' = woi, 



28 - n (2) 



d 



d 



2<5„ P^-P^ 



fly — , m y 

do 



do 



two- level medium is described by the version (©-([7]), 
((12)1. OH of the RMB equations. We believe this state- 
ment to be misleading since the derivation of Eqs. (|12|) . 
(|T5]) contains essential omission. 

To explain why this statement in |l| is incorrect we 
consider the Maxwell equations for the components E' x 
and E' y of the electric field. So, we have 



An evolution of the electromagnetic field of the pulse 
has to obey the Maxwell equations if the semiclassi- 
cal approach is applied. It was claimed in [l[ that the 
transformed components E x and E y satisfy the " Maxwell 
equations" 
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where n is the refractive index of the medium, c is the 
light speed in free space. The quantities P x and P y in 
the right hand sides of Eqs. (© and (© were interpreted 
as the components of the medium polarization and were 
defined for this reason in the following manner: 



P, 



-NTr 



P y = —NTr 




(10) 



(11) 



where N is the density of the medium. 
At last, the reduced equations 



d£ x 
d X ' 
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Rz£y 
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(12) 
(13) 
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(14) 



(15) 



where the components of the medium polarization are 

dH 1 



PL = —NTr 



P' 

v 



-NTr 




(16) 



(17) 



(compare with Eqs. (flQ)). pi)) ). Substitution of ((2j into 
(IT5|) . ((17)) leads to the standard formulas: 

P' x = N Tr (pd x ), Py — NTr (pd y ). 

From Eqs. p"4")) . (TT5)) and the transformation ©, one 
obtains the wave equations on transformed components 
E r , and E v : 
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(18) 
(19) 

where quantities P x and P y are defined by the relation 



where 



P, 



= T 



P' 
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P' 
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(20) 



R, = 



nhc V n ) 



were obtained from ((H)) and @ with the help of the uni- 
directional propagation approximation. 

The misprints were corrected in the formulas presented 
above. We divided the first term in the right hand side 
of pi by two, changed a sign in the definition of fj, y , 
multiplied variable x' by wo/2. Also, Eqs. |8))- pT)) were 
written in the terms of the variables E x , E y instead of 
£ x , £ y , and the multipliers d x , d y are omitted in the 
right hand sides of Eqs. (O, ©■ We found that these 
corrections are necessary for the system (0- ((7)), p2)). 
(JT3J) to be obtained. 

According to Ref. [l| , an evolution of the unidirectional 
two-component electromagnetic pulses in an anisotropic 



A connection between the components P x , Py and P x , 
P y exists also. Taking into account ((3)), we deduce from 
Eqs. (dSJ), CP and (JTB)), p7)) that 
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Define the dimensionless parameter 
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If the condition 



S X Sy 



e = 



(21) 



(22) 



holds, then matrix T is orthogonal: T T — T 1 . It can 
easily be seen from Eqs. (|2"U)) and (|2ip that P. x = P x , 
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P y = Py in this case. Formula ([3]) under this condition 
is nothing but the rotation transformation. 

In the general case (e ^ 0), we have from (|20| . (|2T|) 
that P x P x and P y ^ P y , i.e. the right hand sides of the 
equations ([8]), © and (TTg| . (fl~9f on transformed compo- 
nents -Ez, -Ej, are different. This shows that it is incorrect 
to determine the quantities P x and P y in Eqs. (JSj> , (|9|) by 
means of the formulas (TIU|) and (fTT]) if e 7^ 0. These for- 
mulas are valid in the particular case e = when trans- 
formation © is the rotation transformation. Thus, the 
version ©-©), HU, (p© of the RMB equations can be 
applied only if the condition (f2"2"]l is imposed on the ele- 
ments of the matrices d x and d y . It is wrong to exploit 
this version in the general case. 

An application of the unidirectional propagation ap- 
proximation to Eqs. (|18[) and (TIT))) gives 



7 — Rz£y — ^xSy — e(R z £ x — /j, y S x ), (23) 



7^7 = —Rz£x + H V S X + e(R z £ y - fi x S y ). (24) 



The wrong version ©-©, 0©, 0© of the RMB equa- 
tions contains four parameters /i x , fi y , m 
nected by the relation 



and rriy con- 



1. 



(25) 



It was claimed in [T[ that the Lax pair exists for this 
version. 

The correct version t[23 |) . t[24 |) of the RMB 

equations contains five parameters fi x , fi y , m x , m y and e. 
The relation (gU) is fulfilled also. The system of the RMB 



equations equivalent to the correct version was consid- 
ered in Ref. It was shown that this system possesses 
the Lax pair in the most general anisotropic case when 
all the elements of the matrices d x and d y are arbitrary. 

The system of the RMB equations and its Lax pair 
were written in Q in the terms of physical variables and 
parameters. Having rewritten these systems in the terms 
of the variables t', x' an d £0 £y, we obtain equations 
©-©, ([M]), J2D and their Lax pair: 



dtp 
dx 1 



A(X)ip. 



(26) 
(27) 



Here tp = -0(t', x', A) is a solution of the Lax pair, A is 
the spectral parameter, 2x2 matrices L(X) and A(X) are 



defined as given 
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(28) 



(29) 



E = 



Ex s*E y S 3 5 2 Q 

7r=^ + ^' R= Q + ^a' 



5i = -Vl 



5 3 S Z + S 4 S X + 5 5 S y , 



s 2 = 



s 2 (n x + ml) + fly + rriy + 2sm x m y 
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It can be checked immediately that the overdetermined 
system (|2"6"|) . (j2"7| is the Lax pair of the correct version 
of the RMB equations ©-©, ((Ml), (EU- Indeed, the 
compatibility condition of Eqs. (|26| . (|27|) is 



8L(X) dA(X) 



dx' 



dr> 



[L(X),A(X)}=0. 



(30) 



A substitution of (US]) and into ©I]) yields the system 
of the RMB equations I®, I®. 

The system equivalent to Eqs. ©-0, ((23]), ([24]) is 
obtained by applying the unidirectional propagation ap- 
proximation to the Maxwell equations (TT4"]) . (|T5|) . The 
coefficients of the Lax pair found in [2| for this system 
were expressed directly through the elements of the ma- 
trices oL and d„. 
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